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INTRODUCTION 


This  report  describes  a  finite  difference  computer  program  for 
investigating  the  response  of  multilayered  gun  barrels  subject  to  some  firing 
cycle.  Results,  typical  of  which  the  program  is  capable  of  generating,  are 
presented  for  a  tube  model  which  has  a  TZM  liner  and  a  steel  jacket.  The 
application  is  to  a  large  caliber  weapon,  but  the  program  can  be  used  for 
small  caliber  also.  The  program  was  written  to  coincide  with  a  development 
program  which  is  examining  the  feasibility  of  fabricating  and  firing 
multilayered  gun  tubes.  One  of  the  main  factors  limiting  the  life  of  gun 
tubes  is  the  excessive  wear  and  erosion  which  occurs  especially  at-  the  forcing 
cone  area  of  the  gun  tube.  The  experimental  program,  which  has  shown  success 
for  20  ram  weapons,  is  to  insert  liners  fabricated  from  refractory  materials 
into  the  forcing  cone  area  of  the  gun  tube.  Since  refractory  materials  have 
high  melting  points,  there  is  a  strong  indication  that  they  will  experience 
less  wear  and  thus  Increase  the  life  of  the  weapon. 

Earlier  versions  of  the  computer  program  have  been  used  to  describe  other 
behavior. 1“* 2 3  Preliminary  work  on  the  current  problem  was  presented  in 
reference  3.  That  work  has  been  improved  by  inputting  the  thermo-physical 
properties  as  functions  of  temperature  and  by  allowing  contact  resistance 


IVasilakis,  J.  D. ,  "Temperatures  and  Stresses  Due  to  Quenching  of  Hollow 
Cylinders,"  Transactions  of  the  24th  Conference  of  Army  Mathematicians,  ARO 
Report  79-1,  January  1979. 

2Vasilakis,  J.  D.  and  Chen,  P.  C.  T. ,  "Therrao-Elastic-Plastic  Stresses  on 
Hollow  Cylinders  Due  to  Quenching,"  Transactions  of  the  25th  Conference  of 
Army  Mathematicians,  ARO  Report  80—1,  January  1980. 

3Vasilakis,  J.  D.  ,  "Thermo-Elastic  Plastic  Stresses  in  Multi-Layered 
Cylinders,"  Transactions  of  the  27th  Conference  of  Army  Mathematicians,  ARO 
Report  82-1,  January  1982. 
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between  layers.  The  boundary  conditions  have  also  been  improved  so  they  are 
now  generated  for  the  current  problem  at  hand,  i.e»,  for  specific 
configuration  and  bore  material,  whereas  previously  they  were  empirically 
generated  for  another  system  and  simply  used  in  the  program  as  typical  input. 

The  computer  program  consists  of  three  parts  which  can  be  run  as  a  single 
program  or  as  three  separate  stand  alone  programs.  The  first  is  an  internal 
ballistics  program  which  generates  the  boundary  conditions,  i.e.,  heat 
transfer  coefficients,  pressures,  and  gas  temperatures  as  a  function  of  time 
for  a  single  firing  pulse,  for  input  to  the  next  two  programs.  The  next 
program  section  computes  the  transient  temperatures  due  to  some  firing  cycle 
and  can  be  used  to  show  the  thermal  response  of  the  system  over  several  firing 
cycles  indicating  the  temperature  buildup,  and/or  the  temperatures  can  be  used 
as  input  to  the  third  program  section  for  the  computation  of  stresses.  This 
can  be  done  as  the  thermo-mechanical  program  is  treated  as  uncoupled. 

Results  from  each  of  the  three  sections  using  a  TZM  liner/steel  jacket 
configuration  for  a  105  mm  large  caliber  weapon  system  are  presented. 

PROCEDURE 

The  equations  used  to  describe  the  behavior  are  discussed  below  followed 
by  the  numerical  work. 

Theory 

The  partial  differential  equation  for  describing  the  axisymmetric 
transient  temperature  distribution  in  multilayered  cylinders  is  given  by ,  for 

layer  L,  T  L 

18  3T  T  T  3T 

-  (kL(T)r  —  )  =  cL(T)pL(T)  —  (1) 

r  3r  3r  3t 
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where  r  represents  the  radial  distance,  T  the  temperature,  and  t  the  time. 

The  thermal  conductivity,  specific  heat,  and  density  are  given  by  k,  c,  and  p 

respectively.  These  properties  are  assumed  to  be  functions  of  temperature. 

Axial  effects  are  ignored  in  the  program.  The  geometry  is  shown  in  Figure  1. 

The  initial  condition  is  given  by 

T(r,o)  =  Tq  (2) 

where  T0  would  normally  represent  some  ambient  temperature.  A  temperature 

other  than  ambient,  say  due  to  some  environmental  condition,  could  also  be 

used.  The  boundary  conditions  are  of  the  type 

8T 

k(T)  —  -  h(T-Tg)  =  -g  (3) 

3r 


where  h  is  convection  type  heat  transfer  coefficient,  g  would  represent  some 
heat  input  if  it  existed,  and  Tg  is  the  temperature  of  the  propellant  gases 
when  the  boundary  condition  is  applied  on  the  inside  or  bore  diameter  and  the 
ambient  temperature  when  applied  on  the  outside  surface  of  the  gun  tube. 

The  thermo -physical  properties  are  made  dimensionless  with  regard  to 
their  respective  values  for  steel  at  the  ambient  temperature.  The  tempera¬ 
tures  are  made  dimensionless  with  regard  to  the  maximum  gas  temperature 
achieved  during  the  interior  ballistic  cycle.  The  dimensionless  time  is 
defined  by 

kQt 

t  - - -  (4) 

Pocob 


where  ko,  pQ,  and  cQ  are  the  values  for  steel  as  mentioned  previously. 


the  boundary  conditions, 

A  hb 
h  =  — 
k 

becomes  the  Nusselt  number. 


In 

(5) 
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Also  required  are  continuity  conditions  between  the  concentric  cylinders. 
Normally  one  requires  continuity  of  temperature  and  heat  flux.  However, 
contact  resistance  does  exist  between  surfaces  and  it  was  decided  to  include 
this  effect  here.  The  resistance  results  from  the  true  nature  of  surfaces.^ 
Conduction  occurs  at  the  discrete  points  of  contact  between  the  surfaces  and 
is  therefore  a  function  of  pressure,  surface  conditions,  fluids  in  the  voids, 
etc.  It  is  treated  here  as  a  thin  layer  resisting  the  flow  of  heat.  For  the 
finite  difference  formulation,  there  is  a  jump  in  temperature  at  node  i 

TiL  _  xiL+1  +  AT  (6) 


while  the  heat  flux 


-  =  k^(T) 
A 


Ti-1L  -  TiL 


Ar 


=  hcAT  =  kL+1(T) 


Ti 


L+l  _ 


Ti+1 


L+l 


Ar 


(7) 


remains  constant  as  one  passes  across  the  layer  from  cylinder  L  to  cylinder 
L+l.-  Equations  (6)  and  (7)  allow  the  computation  of  AT  and  TxL+1. 

The  use  of  finite  difference  equations  to  solve  the  thermo-elastic- 
plastic  stress  problem  requires  expressing  the  equilibrium  equation  and  the 
equation  of  compatibility  at  each  node  at  which  the  finite  difference 
equations  are  desired.  The  Prandtl-Reuss  flow  rule  is  used  to  eliminate  the 
incremental  stresses  so  that  what  results  is  a  matrix  for  evaluating  the 
incremental  radial  and  tangential  strains  at  each  node.  The  required 
equations  follow,  written  in  dimensionless  form.  The  problem  is  treated  as 
plane  strain. 


4Fenech,  H.  and  Rohsenow,  W.  M. ,  "Prediction  of  Thermal  Conductance  of 
Metallic  Surfaces  in  Contact,"  Journal  of  Heat  Transfer,  February  1963,  pp. 
15-24. 
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The  equation  of  equilibrium  is  written 


where 


3or  ar-ae 

- + - -  0 

3r  r 


is  the  dimensionless  radial  stress 


(8) 


ae 

00  (=  — )  is  the  dimensionless  tangential  stress 


and  oQ  is  the  yield  stress  in  tension,  and  the  compatibility  equation 

8ee  ee-er 

- + - -  0  (9) 

8r  r 


where 


£0 

.  £e(=  E  — )  is  dimensionless  tangential  strain 

°o 


er(=*  E  — )  is  dimensionless  radial  strain 


and  o0/E  is  yield  strain  in  tension  when  E  is  Young’s  Modulus.  The 
compressibility  of  the  material  is  expressed  by 

o 

e  -  aT  +  — 

3K 

1 

e  ■  -  (er+ee)  is  mean  strain 
1 

a  *  -  (or+O0+az)  ^an  stress 


(10) 


K 

K(=  — )  is  dimensionless  bulk  modulus 
°o 

a(“  aT^)  is  dimensionless  coefficient  of  thermal  expansion 
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and 

ez  ~  0  for  plane  strain 

Traction  free  boundary  conditions  are  always  used  in  the  outside  radius 
and  on  the  bore  when  only  thermal  stresses  are  required.  When  mechanical 
loads  are  desired,  the  pressure  pulse  is  applied  to  the  bore. 

It  was  desirable  to  write  the  finite  difference  equations  in  terras  of 
strain  alone;  hence,  the  stresses  in  the  equations  of  equilibrium  had  to  be 
expressed  in  terras  of  the  strains.  This  was  accomplished  by  modifying  a 
plastic  stress-strain  matrix5  which  was  derived  by  inverting  the  Prandtl-Reuss 
equations.  The  inverted  Prandtl-Reuss  equation  is 

„  EodT 

{da}  =  [Dp](de} - {1}  (11) 

(l-2v)a0 

where  the  stress  vector  is  {da}  =  {dar,  dag,  doz}T,  the  strain  vector  {de}  = 
{der,  deg,  0}^,  and  {1}  represents  a  unit  vector.  The  plastic  stress-strain 
matrix  [Dp]  is  given  by 


1— v  °r ' 2 

l-2v  S 


[DP]  =  --- 
1+v 


v  ar'ag' 

l-2v  S 


1-v  00' 2 

l-2v  S 


v  °r ' °z '  v  og' az ' 

l-2v  S  l-2v  S 


V 
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1-v  oz’2 

l-2v  S 


(12) 


^Yamada,  Y.  ,  Yoshimura,  N. ,  and  Sakuri,  T. ,  "Plastic  Stress-Strain  Matrix  and 
Its  Application  For  the  Solution  of  Elastic-Plastic  Problems  by  the  Finite 
Element  Method,"  International  Journal  of  Mechanical  Sciences,  1968,  Vol.  10, 
pp.  343-354. 
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The  primed  stresses  are  deviator ic  stresses. 


ai'  =  °i  “  “  I  01  1  =  r,0,z 


(13) 


At  each  node  during  a  computation,  the  von  Mises1  yield  criterion 

-  l(or-O0)2  +  (a0-ar)2  +  (az-ar)2]  =1  (14) 

2 

is  checked  to  see  if  plastic  deformation  has  progressed  to  that  node.  If  not, 
the  stresses  remain  elastic  and  can  still  be  computed  using  Eq.  (12)  by 
setting  the  deviatoric  stresses  equal  to  zero.  The  matrix  [Dp]  then  becomes 
the  same  matrix  as  would  exist  if  linear  elastic  behavior  had  been-  assumed. 
The  quantity  S  is  given  by 


where 


S 


(1 


-3  3 

o  =  -  Oij'aij  =  -  (a 


t  2 


r 


00 


» 2 


>Z2) 


(15) 


(16) 


is  the  equivalent  stress  and 


do 

H'  - - 

den 


(17) 


is  the  slope  of  the  equivalent  stress/equivalent  plastic  strain  curve  and  is  a 
measure  of  hardening.  The  increment  in  equivalent  plastic  strain  is  given  by 


dep  =  -  deijPdeijP 


08) 
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Boundary  Conditions 


The  boundary  conditions  that  are  used  as  input  for  the  calculation  of 
temperatures^  and  of  stresses  are  generated  using  a  computer  program  based  on 
reference  6.  In  that  paper,  the  burning  of  a  specific  propellant  for  the 
purpose  of  firing  a  projectile  from  a  gun  tube  is  modeled.  The  equations  used 
are  based  on  Corner’s  work7  and  represent  a  first  order  interior  ballistics 
solution.  Lagrange’s  approximation  is  assumed,  i.e.,  the  velocity  of  the  gas 
at  any  instant  increases  linearly  with  distance  along  the  bore  from  zero  at  the 
breech  to  the  full  shot  velocity  at  the  base  of  the  projectile.  Exponential 
decay  is  assumed  during  the  blowdown  cycle,  i.e.,  after  the  projectile  has  left 
the  gun  tube.  Based  upon  the  rate  equations  and  heat  balances  involved,  the 
heat  transfer  coefficients,  the  pressure  pulse,  and  the  gas  temperature  can  be 
found  as  a  function  of  time  during  the  firing  cycle.  The  quantity  of  heat  that 
goes  into  the  heating  of  the  gun  tube  can  be  computed.  The  bore  surface  can  be 
specified  to  be  the  liner  material.  Figure  2  shows  the  output  of  this  program 
for  a  gun  tube  with  steel  at  the  bore  surface. 

Numerical  Procedure 

The  Crank-Nicolson  representation  for  finite  differences  of  the  partial 
differential  equation  governing  the  temperatures  in  time  is* 


^-Vasilakis,  J.  D.  ,  "Temperatures  and  Stresses  Due  to  Quenching  of  Hollow 
Cylinders,*’  Transactions  of  the  24th  Conference  of  Army  Mathematicians,  ARO 
Report  79-1,  January  1979. 

^Vottis,  P.  M. ,  "Digital  Computer  Simulation  of  the  Interior  Ballistic  Process 
in  Guns,"  WVT-6615,  Watervliet  Arsenal,  Watervliet,  NY,  October  1966. 

7 Corner,  J.  ,  Theory  of  Interior  Ballistics  of  Guns,  John  Wiley  and  Sons,  Inc., 
New  York,  1950. 
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[a+iAr)ki+i/2,n+l/2lTi+l ,n+l  +  t“(a+iAr)ki+l/ 2 ,n+l/ 2 

2Ar2 

-(a+(i-l)Ar)ki_1/2)n+l/2pi,n+l/2(““)(a+(i-1/2)Ar)Ti,n+l 
+  [(a+(i-l)Ar)ki_1/2.n+l/2]Ti-l,n+l  “  [-(a+iAr)k1+i/ 2 ,n+l/ 2l Ti+1 ,a  + 


+  [(a+iAr)ki+i/2,n+l/2  +  (a+( i-1) Ar)ki_i/ 2 ,n+l/ 2  “ 

2Ar2 

Ci,n+l/2pi,n+l/2(----)(a+(i-1/2>Ar)lTi,n  +  l-(a+(i-l)Ar)ki_1/2,n+l/2JTi-l ,n 

(19) 

The  equation  is  solved  twice: 

1.  At  n+1/2  step,  allowing  k,p,c  etc.  to  take  on  the  values  at  t  =  n 

step. 

2.  The  new  temperatures  are  then  used  to  evaluate  k,c,p  at  n+1/2  step  and 
the  set  of  equations  re-evaulated  for  the  temperatures  at  the  n+1  step. 

The  computed  temperature  distributions  at  each  full-time  step  are  saved  on 
disk  and  eventually  called  in  when  required  by  the  stress  program. 

The  finite  difference  equations  are  within  any  layer. 

Compatibility : 

-rjAce^  +  (2ri-ri-i)Aee^  -  (ri-ri-i)Aer^  = 


-r1(ee^-ee^  -  (ri-ri-i) ( ee^-ee^ ) 


(20) 


Equilibrium: 


-r-iAar  -  (ri-ri-i)Aae  +  (2ri-ri-i)Aar 
i-1  i  1 

-ri(°r  “°r  , )  “  (ri“ti-i)(ar  -00  ) 

i  i-1  i  i 


(21) 
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Substituting  the  Prandtl-Reuss  equations  into  that  of  equilibrium 

-r ^D(r , 0 )Ae -riD(r ,r ) Aer^  ^  +[-(ri-ri-i)D( 0 , 0)  +  (  2ri-r-L-i)D(r  ,  0)  ]  Agq^ 

+  [-(ri-ri-^DCOjr)  +  (2ri-ri-i)D(r ,r) ] Aer^ 

Ecx 

t±[ot  -ar  ]  +  (ri~ri-i)(ae^-ar^)  +  ri  ~~  [ AT^-ATi-^]  (22) 

Equations  (20)  and  (21)  are  in  backward  difference  equations.  The  actual 
computations  are  performed  by  averaging  backward  and  forward  difference 
schemes.  At  the  interface  between  cylinders,  continuity  of  the  radial  stress 
and  radial  displacement  is  specified  and  on  the  boundary,  i  =  1, 

EaATi 

D(r,0)Ae0  +  D(r,r)Aer  = - -  Ap^  (23) 

1  1  l-2v 

where  Apj_  represents  a  pressure  increment  at  the  bore  or  inside  diameter. 

The  solution  procedure  for  the  transient  temperature  problem  is  as 
follows.  The  temperature  problem  is  solved,  and  the  temperature  distributions 
at  their  computation  times  are  stored  on  disk.  These  distributions  are  called 
into  the  thermo-elastic-plastic  stress  program  one  at  a  time.  The  correspond¬ 
ing  thermal  stresses  are  calculated  and  each  node  checked  to  see  if  the  yield 
criterion  is  satisfied.  If  not,  the  problem  is  still  assumed  to  be  elastic,  a 
new  temperature  distribution  is  called  in,  and  new  stress  increments  calcu¬ 
lated.  The  stresses  are  updated,  and  the  yield  criterion  checked  again.  When 
the  stresses  at  a  point  are  found  to  satisfy  the  yield  criterion,  the  node  is 
identified,  and  the  stress  increments  at  that  node  from  the  next  set  of  temper¬ 
atures  are  computed  using  the  Prandtl-Reuss  equation  or  [Dp]  matrix  identified 
earlier.  This  procedure  is  continued  with  new  sets  of  temperature  called  in 
and  with  the  tracking  of  the  elastic-plastic  boundary  with  time. 
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The  mechanical  properties  are  evaluated  at  the  existing  temperatures. 
However,  the  yield  stress  has  not  yet  been  incorporated  as  a  function  of 
temperature  in  the  program. 

Mechanical  and  Thermo-Physical  Properties 

The  properties  used  in  the  calculations  were  found  in  reference  8.  The 
nominal  values  are  given  in  Table  I.  It  is  always  one  of  the  more  difficult 
tasks  to  find  properties  as  functions  of  temperature.  The  steel  properties 
used  were  those  of  4340  and  4150.  “Gun  Steel"  is  typically  4340  or  a 
modification  thereof.  The  thermo-physical  properties  for  TZM  used  were  those 
for  molybdenum  itself  since  they  were  readily  available  as  functions  of 
temperature  and  the  same  properties  for  TZM,  only  given  at  one  or  two  specified 
temperatures  tended  to  fall  on  or  near  the  same  property  for  molybdenum. 


TABLE  I.  ROOM  TEMPERATURE  PROPERTIES 


k 

c 

P 

a 

E 

°yp 

V 

BTU/#in°F 

BTU/#°F 

#/ in3 

in/in°F 

Psi 

Psi 

Steel 

5.01xl0-lt 

.105 

.289 

6.2xl0-6 

30xl06 

160xl03 

.3 

TZM 

1.87xl0-4 

.06 

.369 

3.0xl06 

45xl06 

130x10 3 

.314 

® Aerospace  Structural  Metals  Handbook,  AFML-TR-68-1 15 . 
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RESULTS 


The  interior  ballistic  code  was  first  run  to  set  up  the  input  data  (heat 
transfer  for  coefficients  and  gas  temperatures  during  firing  cycle)  for  the 
position  of  the  program  which  computes  the  transient  temperature  distribution 
and  pressure  time  curve  for  the  mechanical  load  contribution  to  the  stress  part 
of  the  program.  Two  data  sets  were  established,  one  for  steel  at  the  bore  and 
one  for  TZM  at  the  bore.  Figure  2  shows  the  results  for  steel.  The  interior 
ballistic  code  computes  these  results  at  several  stations  along  the  tube,  but 
only  the  section  at  which  the  maximum  pressure  was  generated  was  considered  at 
this  time.  The  temperature  portion  of  the  program  serves  two  separate  pur¬ 
poses.  The  program  can  be  run  over  several  firing  pulses  oased  on  some 
specific  firing  cycle.  This  will  show  the  buildup  in  temperature  during  the 
firing.  It  can  also  be  run  to  provide  input  data  in  the  form  of  temperature 
distributions  throughout  the  wall  of  the  tube  for  specific  times  during  the 
firing  cycle.  This  data  set  is  then  used  in  the  stress  program  for  the  compu¬ 
tation  of  thermal  stresses  or  thermo -mechanical  stresses  when  the  pressure-time 
curve  is  also  applied.  There  were  four  types  of  material  problems  considered, 
a  single  (monobloc)  steel  tube  with  constant  properties,  monobloc  steel  tube 
with  temperature  dependent  properties,  a  TZM  liner/ steel  jacket  with  constant 
properties  in  each  cylinder,  and  a  TZM  liner/steel  jacket  with  temperature 
dependent  properties.  Results  from  some  of  these  cases  are  presented  below. 

Figure  3  shows  the  change  of  the  bore  temperature  with  time  over  four 
firing  cycles.  The  configuration  is  the  TZM  liner/steel  jacket  with  tempera¬ 
ture  dependent  properties.  The  firing  cycle  depicted  represents  a  projectile 
being  fired  at  the  rate  of  four  rounds  per  minute.  The  temperature  buildup  at 
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the  bore  can  easily  be  seen.  Figure  4  shows  the  temperature  response  at  the 
bore  for  a  monobloc  steel  tube  with  temperature  dependent  properties.  This  is 
shown  on  an  expanded  scale  and  represents  the  thermal  response  used  on  one  of 
the  data  sets  for  the  stresses.  When  these  data  sets  for  stresses  are 
established,  a  finer  time  increment  is  used  over  that  which  simply  computes  the 
thermal  response  for  temperature  buildup.  The  finer  time  increment  decreases 
the  temperature  difference,  AT,  from  time  step  to  time  step  which  is  used  in 
the  stress  program.  This  helps  in  approaching  yield  in  the  stress  program  even 
though  the  temperature  difference  is  further  divided  near  yielding  and  after 
yielding  has  begun.  When  only  the  temperatures  themselves  are  desired,  it  was 
previously  shown 3  that  larger  time  increments  can  be  used.  As  the  finite 
difference  program  for  temperatures  is  implicit,  the  time  increment  between 
pulses  is  increased  at  a  rapid  rate  until  the  next  firing  pulse  comes  along. 

It  took  approximately  90  time  steps  to  complete  the  temperature  response  to  the 
four  cycles. 

When  investigating  the  effect  of  contact  resistance,  it  becomes  obvious 
why  many  papers  ignore  it.  The  difficulty  is  not  with  the  computation,  but 
rather  the  uncertainty  of  the  physical  constant  to  use  in  the  evaluation.  The 
property  h  described  earlier,  is  treated  as  a  constant  here  but  in  reality 
would  be  a  function  of  pressure,  temperature,  the  roughness  of  the  surfaces  in 
contact,  the  hardness  of  the  materials  involved,  etc.  One  would  actually  need 
the  true  area  of  contact  as  opposed  to  the  apparent  area  and  how  this  changed 

^Vasilakis,  J,.  D.  ,  "Thermo-Elastic-Plastic  Stresses  in  Multi-Layered 
Cylinders,"  Transactions  of  the  27th  Conference  of  Army  Mathematicians,  ARO 
Report  82-1,  January  1982. 


13 


in  time.  Table  II  shows  the  effect  of  varying  h  on  the  bore  temperature  for 
five  firing  pulses.  The  h  is  dimensionless.  The  table  shows  that  when  h  = 
1000,  the  system  is  equivalent  to  having  no  resistance  to  heat  flow  and  when 
h  =  .00001,  the  effect  is  equivalent  to  zero  heat  flow  at  the  interface.  As 
can  be  seen  by  converting  some  of  the  values  from  reference  4,  one  can  get 
coefficients  which  result  in  measurable  effects.  To  better  show  the  effect, 
the  resulting  temperature  distribution  as  a  function  of  radius  as  shown  in 
Figure  5  for  h  =  1.  The  second,  third,  and  fourth  pulses  are  displaced  from 
the  first  for  clarity  of  viewing.  The  bore  temperature  increases  substantially 
compared  to  zero  contact  resistance  (as  can  also  be  seen  from  Table  II).  The 
temperature  jump  also  increases  substantially.  The  line  indicating  the  jump 
should  be  a  vertical  line  from  the  point  on  the  inner  cylinder.  The  fact  that 
it  is  not  is  due  to  the  plotting  routine.  A  final  remark  on  this  section  would 
be  that  with  the  uncertainty  in  computing  or  experimentally  determining  an 
actual  h  for  a  system,  Table  II  shows  little  difference  between  an  h  =  1000  and 
h  =*  100  as  a  threshold  for  zero  contact  resistance  and  again  little  difference 
between  h  =*  •  1  and  h  =  .00001  for  an  insulating  barrier. 


^Fenech,  H.  and  Rohsenow,  W.  M. ,  "Prediction  of  Thermal  Conductance  of  Metallic 
Surfaces  in  Contact,"  Journal  of  Heat  Transfer,  February  1963,  pp.  15-24. 
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TABLE  II .  EFFECT  OF  H  ON  BORE  TEMPERATURE 


1000 

100 

1 

.1 

.00001 

Pulse 

i 

0 

0 

0 

0 

0 

Pulse 

2 

.00352 

.00356 

.00837 

.00999 

.01021 

Temp 

Pulse 

3 

.00580 

.00582 

.01493 

.01952 

.02020 

Pulse 

4 

.00758 

.00756 

.02013 

.02861 

.02998 

Pulse 

5 

.00904 

.00901 

.02424 

.03728 

.03596 

The  remaining  results  relate  to  the  output  from  the  stress  portion  of  the 
program •  Figure  6  shows  the  effect  of  including  temperature  dependence  of 
material  properties  in  the  analysis.  The  variation  of  radial  stress  and 
tangential  stress  across  the  wall  thickness  is  shown.  There  is  little  effect 
except  at  the  peak  tangential  stress  at  the  bore  which  is  about  three  to  four 
percent  higher  for  a  tube  with  temperature  dependent  material  properties.  This 
represents  a  stress  difference  of  about  5000  psi. 

The  stress  results  being  shown  are  only  the  response  to  the  first  firing 
pulse.  Hence,  the  stresses  are  mainly  due  to  a  steep  thermal  gradient  at  the 
bore  of  the  tube.  If  inelastic  response  does  not  occur,  the  stress  response 
would  not  change  significantly  from  pulse  to  pulse  as  the  mechanical  load 
vanishes  and  there  remains  only  a  slight  thermal  gradient  throughout  the  tube 
wall.  Depending  on  the  configuration  of  the  system,  the  firing  rate,  etc., 
the  rise  in  bore  temperature  before  the  next  round  is  fired  is  small,  ~  3°F  to 
9°F . 
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Figure  7  shows  the  variation  in  tangential  stress  versus  radius  for  the 
multi-cylinder  configuration,  TZM/steel,  with  the  temperature  dependent  proper¬ 
ties  at  the  time  when  the  internal  pressure  peaks.  The  mechanical  loading  due 
to  the  pressure— t ime  curve  alone  is  shown  as  is  the  thermal  stress  distribution 
due  to  the  temperature  distribution.  The  combined  effects  are  also  shown. 

These  combined  effects  are  not  arrived  at  by  adding  the  separate  ones  but  are 
recomputed.  This  is  important,  especially  in  the  case  shown,  because  the 
mechanical  load  alone  was  sufficient  to  cause  the  inelastic  deformation  at  the 
bore.  While  this  inelastic  zone  was  concentrated  at  the  bore  and  there  was 
little  depth  of  plastic  zone  penetration  into  the  wall,  it  was  still  incorrect 
to  linearly  superimpose  solutions  by  adding  the  separate  stress  distributions. 
Under  the  combined  loads,  the  solution  remained  elastic  throughout.  Figure  3 
shows  the  results  from  the  same  problem  at  peak  thermal  load. 

SUMMARY 

The  computer  program  described  is  capable  of  predicting  the  thermal 
response  and  the  thermo-elastic-plastic  response  of  liner/ jacket  gun  tube 
designs  with  temperature  dependent  thermo— physical  properties.  Realistic  input 
to  this  problem  is  generated  using  an  interior  ballistics  code.  While  only  a 
two  cylinder  system  is  described,  allowance  is  made  for  up  to  five  cylinders. 

Improvements  can  always  be  made  and  in  this  case,  the  following  could  be 
included.  By  using  a  variable  space  increment  in  the  thermal  program,  the 
effect  of  a  deposited  thin  layer  on  the  bore  surface  can  be  investigated. 
Initial  stresses  due  to  interference  should  be  Incorporated  as  should  a 
temperature  dependent  yield  stress. 
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Figure  1.  Typical  Multi-Layered  Geometry 
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Figure  2, 
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Figure  3.  Bore  Temperature  Response  over  Several 
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Figure  6.  STRESS  v'S  RADIUS 
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